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ABSTRACT 

We study the distribution of dark matter in dwarf spheroidal galaxies by modelling 
the moments of their line-of-sight velocity distributions. We discuss different dark 
matter density profiles, both cuspy and possessing flat density cores. The predictions 
are made in the framework of standard dynamical theory of two-component (stars 
and dark matter) spherical systems with different velocity distributions. We compare 
the predicted velocity dispersion profiles to observations in the case of Fornax and 
Draco dwarfs. For isotropic models the dark haloes with cores are found to fit the 
data better than those with cusps. Anisotropic models are studied by fitting two 
parameters, dark mass and velocity anisotropy, to the data. In this case all profiles 
yield good fits but the steeper the cusp of the profile, the more tangential is the 
velocity distribution required to fit the data. To resolve this well-known degeneracy 
of density profile versus velocity anisotropy we obtain predictions for the kurtosis 
of the line-of-sight velocity distribution for models found to provide best fits to the 
velocity dispersion profiles. It turns out that profiles with cores typically yield higher 
values of kurtosis which decrease more steeply with distance than the cuspy profiles, 
which will allow to discriminate between the profiles once the kurtosis measurements 
become available. We also show that with present quality of the data the alternative 
explanation of velocity dispersions in terms of Modified Newtonian Dynamics cannot 
yet be ruled out. 

Key words: methods: analytical - galaxies: dwarf - galaxies: fundamental parameters 
- galaxies: kinematics and dynamics - cosmology: dark matter 


1 INTRODUCTION 

Dark matter is now generally believed to be the dominat¬ 
ing component of the matter content of the Universe which 
has played an important role in the formation of struc¬ 
ture. Significant effort has gone into establishing its prop¬ 
erties and mechanisms of evolution. The theoretical results 
which got most attention in the recent years concerned the 
shape of density profiles of dark matter haloes. Those were 
mainly obtained by the means of TV-body simulations with 
the most popular being the so-called universal profile advo¬ 
cated by Navarro, Frenk and White (1997, hereafter NFW). 
NFW found (and others confirmed) that in a large range 
of masses the density profiles of cold dark matter haloes 
forming in various cosmologies can be fitted with a simple 
formula with only one fitting parameter. This density pro¬ 
file steepens from r -1 near the centre of the halo to r~ 3 at 
large distances. More recent simulations with higher resolu¬ 
tion produce even steeper density profiles, with inner slopes 
r -3 / 2 (Fukushige & Makino 1997; Moore et al. 1998; see also 
Jing & Suto 2000). 

Analytical calculations also predict density profiles with 


steep inner cusps. The NFW profile is claimed to be repro¬ 
duced in studies taking into account the merging mecha¬ 
nism (Lacey & Cole 1993) in the halo formation scenario 
(NFW, Avila-Reese, Firmani & Hernandez 1998). On the 
other hand, using an extension of the spherical infall model 
Lokas (2000) and Lokas & Hoffman (2000) have shown that 
if only radial motions are allowed the inner profiles have to 
be steeper than r~ 2 (in agreement with considerations based 
on distribution functions, see Richstone & Tremaine 1984). 
Although addition of angular momentum is known to make 
profiles shallower (White & Zaritsky 1992; Hiotelis 2002), 
its origin is not yet clear enough to make exact predictions. 

The expected steep inner profiles are now being tested 
by comparison with observations of density profiles of galax¬ 
ies and galaxy clusters. Some recent studies of clusters (e.g. 
van der Marel et al. 2000) claim good agreement between 
cluster observations and the NFW mass density profile, 
while other based on lensing find the dark matter to exhibit 
a smooth core (Tyson, Kochanski & dcll’Antonio 1998, see 
also Williams, Navarro & Bartelmann 1999). In the case of 
galaxies, however, the situation is even more troubling. Flo¬ 
res & Primack (1994) and Moore (1994) were the first to 
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indicate that the NFW profile is incompatible with the ro¬ 
tation curves of spiral galaxies. More and more evidence is 
being presented now, which shows that spiral galaxies pos¬ 
sess a dark matter core rather than a cusp in the inner parts 
(Burkert 1995; Salucci & Burkert 2000; Borriello & Salucci 
2001). The evidence is especially convincing in the case of 
LSB galaxies which are supposed to be dominated by dark 
matter (de Blok et al. 2001; de Blok, McGaugh & Rubin 
2001). Some mechanisms to explain the formation of such 
core have also been proposed (Hannestad 1999; Kaplinghat, 
Knox & Turner 2000; El-Zant, Shlosman & Hoffman 2001; 
Ravindranath, Ho & Filippenko 2001). However, as empha¬ 
sized by van den Bosch & Swaters (2001) even good qual¬ 
ity rotation curves cannot always discriminate between con¬ 
stant density cores and cusps because of the effects of beam 
smearing, the uncertainties in the stellar mass-to-light ratio 
and the limited spatial sampling of the halo density distri¬ 
bution. 

Dwarf spheroidal (dSph) galaxies provide a unique test¬ 
ing tool for the presence and distribution of dark matter be¬ 
cause due to their large velocity dispersions they are believed 
to be dominated by this component (for a review see Mateo 
1997). Measurements of central velocity dispersion allowed 
many authors to estimate mass-to-light ratios in dwarfs and 
find it to be much larger than stellar values (see Mateo 1998a 
and references therein). These calculations were made with 
the simplifying assumption of isotropic velocity distribution 
and in most cases of mass following light. With the mea¬ 
sured velocity dispersion profiles becoming available now, 
this assumption can in fact be relaxed and much more in¬ 
formation can be extracted from the data. However, as in all 
pressure supported systems, when trying to model the ve¬ 
locity dispersion dependence on distance we are faced with 
the well-known degeneracy of density profile versus velocity 
anisotropy (Binney & Tremaine 1987). 

The purpose of this work is to study moments of the 
line-of-sight velocity distribution following from different 
dark matter density profiles and different assumptions con¬ 
cerning velocity anisotropy. We focus on modelling the ve¬ 
locity dispersion profiles and extend here the earlier work 
(Lokas 2001), which discussed only NFW profile, by consid¬ 
ering both cuspy profiles and those possessing flat cores, in 
order to decide which of them (and under what conditions) 
are favoured by observational data in the case of dSph galax¬ 
ies. In order to discriminate between different models we 
also obtain predictions for the fourth-order velocity moment. 
Our work is complementary to the approach of Wilkinson 
et al. (2002) and Kleyna et al. (2001, 2002) who discussed 
dark matter distributions differing at large distances from 
the centre. 

The paper is organized as follows. In Section 2 we briefly 
present the method of calculation of the line-of-sight velocity 
moments in models with different density profiles and veloc¬ 
ity anisotropy. Section 3 describes the assumptions about 
the matter content of dSph galaxies. Results for the veloc¬ 
ity dispersion profiles of Fornax and Draco dwarfs for both 
isotropic and anisotropic velocity distribution are given in 
Section 4, where we also present predictions for the fourth 
moment, the kurtosis. Section 5 provides comments on the 
possible alternative explanation of velocity dispersions in 
terms of Modified Newtonian Dynamics (MOND) and the 
discussion follows in Section 6. 


2 MOMENTS OF THE VELOCITY 
DISTRIBUTION 

The dynamics of stars in a gravitational potential <f>(r) is 
completely described by the distribution function /( r, v) 
governed by the Boltzmann equation. It is often advanta¬ 
geous however to work not with / but with its velocity mo¬ 
ments defined as 

vv l r v J e v*= J v l r v 3 e vlf(r,v) d 3 v, (1) 


where v(r ) is a 3D density of stars. In the following we will 
assume that the system is spherically symmetric and that 
there are no net streaming motions (e.g. no rotation). Then 
all the odd velocity moments vanish. 

At second order the two distinct moments are Vr and 
Vg = which we will denote hereafter by a 2 and oj respec¬ 
tively. They are related by the lowest order Jeans equation 
derived from the Boltzmann equation (Binney & Tremaine 
1987) 


d / 2\ . 2/3 2 . d4? 

— (zzcr r ) -I - + V— = 0, 

dr r dr 


( 2 ) 


where 


0 = 1 - 


(j?(r) 


(3) 


is a measure of the anisotropy in the velocity distribution 
and 4? is the gravitational potential. 

The most convenient approach to solving equation (|i|) 
for of is to make an assumption about /3. The values of this 
parameter have been studied both observationally and via 
V-body simulations but only in systems very different from 
dSph galaxies. For dark matter haloes Thomas et al. (1998) 
find that, in a variety of cosmological models, the ratio ae /cr r 
is not far from unity and decreases slowly with distance from 
the centre to reach ~ 0.8 at the virial radius. Observations 
of elliptical galaxies are also consistent with the orbits being 
isotropic near the centre and somewhat radially anisotropic 
farther away, although cases with tangential anisotropy are 
also observed (Gerhard et al. 2001). In our considerations 
the anisotropy parameter describes the properties of the ob¬ 
served motion of stars which are only a tracer population 
of the underlying dark matter potential and therefore we 
cannot assume their anisotropy is similar to that of dark 
haloes or elliptical galaxies. In the absence of theory of their 
formation, which could provide some hints, or direct mea¬ 
surements of velocity anisotropy in dSph galaxies we have 
to consider different values of j3. 

Traditionally, /3 has been modelled in a way proposed by 
Osipkov (1979) and Merritt (1985) where (3om = r 2 /(r 2 + 
r a ), with r a being the anisotropy radius determining the 
transition from isotropic orbits inside to radial orbits out¬ 
side. This model covers a wide range of possibilities from 
isotropy (r a —> oo,/3 = 0) to radial orbits (r a —> 0,/3 = 1). 
Another possibility is that of tangential anisotropy with the 
fiducial case of circular orbits when /3 —> — oo. For our pur¬ 
poses here it is most convenient to cover all possibilities from 
radial to isotropic and circular orbits with a simple model 
of /3 =const and —oo < (3 < 1. 

The solution of the Jeans equation (|i|) with the bound¬ 
ary condition <r r —> 0 at r —> oo for /3=const is 
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vcr 2 ([3 = const) = r 2/3 f r 2,3 u^- dr. (4) 

J r dr 

From the observational point of view, an interesting, mea¬ 
surable quantity is the line-of-sight velocity dispersion ob¬ 
tained from the 3D velocity dispersion by integrating along 
the line of sight (Binney & Mamon 1982) 


ofos {R) 



\fr 2 — R 2 


dr 


(5) 


where I(R) is the surface brightness. 

Introducing result (Q) into equation (0) and inverting 
the order of integration the calculations of <ti os can be re¬ 
duced to one-dimensional numerical integration of a formula 
involving special functions for arbitrary /3 =const. 

It has been established that by studying cri os (f?.) alone 
we cannot uniquely determine the properties of a stellar sys¬ 
tem. In fact systems with different densities and velocity 
anisotropies can produce identical a\ os (R) profiles (see e.g. 
Merrifield & Kent 1990; Merritt 1987). As will be shown 
below, the same kind of degeneracy also plagues the mod¬ 
elling of dSph galaxies. It is therefore interesting to consider 
higher-order moments of the velocity distribution hoping 
that the models with the same <Ji os (R) will differ in their 
projected higher-order moments. 

In the case of fourth-order moments we have three dis¬ 
tinct components Vr, Vg = v | and v 2 Vg = v 2 v^ related by 
two higher order Jeans equations (Merrifield & Kent 1990). 
In order to calculate the moments we need additional infor¬ 
mation about the distribution function. According to Jeans 
theorem the anisotropic distribution functions can be con¬ 
structed as functions of the integrals of motion, the energy 
E = —4? — v 2 /2 and angular momentum L = (Vg + w^) 1 ^ 2 r. 
In the isotropic case the distribution function is a function 
of energy alone and can be determined uniquely from the 
underlying mass distribution. However, as discussed by De- 
jonghe (1987), in the anisotropic case, in which we will be 
most interested here, the function is not uniquely deter¬ 
mined by the assumed mass density, even if information on 
velocity dispersions is available. 

Therefore we will restrict ourselves here to functions 
which can be constructed from the energy-dependent distri¬ 
bution function by multiplying it by some function of angu¬ 
lar momentum. We will adopt the form 

f(E,L) = ME)L~ 213 (6) 


with j3 =const, a generalization first considered by Henon 
(1973). The velocity moments of such a function can be ex¬ 
plored without explicit calculation of fo{E). For example, 
using equation ([3]) one can show that the function of the 
form (j^) has the desired property of ag /of = 1-/3. Similar 
calculation proves that the three fourth-order moments are 
related by 

Vri>g = i(l-/3)u/ (7) 

V e = ~(l-/3)(2 —/?K. (8) 


Then the two Jeans equations for the fourth-order moments 
reduce to one of the form 


J-(uVr) + —VVr + 3t/of = 0, 
dr r dr 


(9) 
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similar to equation (p|). The solution, in analogy to equation 

(0), is 

Wr{d = const) = 3r -2 ^ / r 2l3 ua 2 (r)^- dr. (10) 

Jr dr 

By projection we obtain the line-of-sight fourth moment 


.(*) 




2 

+ 2(1-0) 


R 2 (r 2 — R 2 ) ( (2 — 0)(1 — 0) R 4 

X -5-1 ^- T 


Vr 2 - R 2 


dr 


( 11 ) 


where we used relations (0m|)- Introducing equations ( |Io| ) 
and (1) into ( pd|) and inverting the order of integration the 
calculation of vf os (R) can be reduced to two-dimensional nu¬ 
merical integration. A useful way to express the fourth pro¬ 
jected moment is to scale it with af os in order to obtain the 
kurtosis 


^los (R) 


< S (R) 

< S (R) 


which will be discussed in Subsection 4.3. 


( 12 ) 


3 MATTER CONTENT 

We will now discuss different mass distributions contributed 
by stars and dark matter to the gravitational potential <t> in 
equation (|i:). 


3.1 Stars 


The distribution of stars is modelled in the same way as in 
Lokas (2001) i.e. by the Sersic profile (Sersic 1968, see also 
Ciotti 1991) 

I(R) = I 0 exp[-(i?/R s ) 1/m ], (13) 


where Io is the central surface brightness and Rs is the char¬ 
acteristic projected radius of the Sersic profile. The best¬ 
fitting Sersic parameter m has been found to vary in the 
range 1 < m < 10 (Caon, Capaccioli & D’Onofrio 1993) for 
different elliptical galaxies, however for dSph systems m = 1 
is usually used, although in some cases other values of m are 
found to provide better fits (e.g. Caldwell 1999). 

The 3D luminosity density v(r) is obtained from I(R) 
by deprojection 


u{r) = 


i r dI — 

7T Jr d R VR 2 - r 2 


(14) 


In the case of m = 1 we get u(r, m = 1) = 

IoKo(r/Rs)/(nRs), where Ko(x) is the modified Bessel 
function of the second kind. For other values of m in the 
range 1/2 < m < 10 an excellent approximation for v(r) is 
provided by (Lima Neto, Gerbal & Marquez 1999) 


/ r \~P 

/ r \ !/ m 

U0 (-R S ) 6XP 

~\R^J 

7 0 r(2m) 


2R s r[(3 -p)m\ 



p = 1.0 - 0.6097/m + 0.05463/m 2 . 


(15) 


The mass distribution of stars following from 


© 


is 
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M,{r) = TLtot 


7[(3~ p)m, (r/R s ) 1/m ] 
r[(3 -p)m] 


(16) 


where T=const is the mass-to-light ratio for stars, Ltot is the 
total luminosity of the galaxy and 7 ( 0 , x) = /o’e-V - 1 di 
is the incomplete gamma function. 


3.2 Dark matter 


We explore here different density distributions of dark mat¬ 
ter which can be described by the following general formula 

PM = (r/r . )n( ?Tr / r. ) a- < 17 > 

where p c h ar is a constant characteristic density. As is im¬ 
mediately clear from equation ([u|), the inner slope of the 
profile is r~ a while the outer one is always r^ 3 . We choose 
the parameter a in the range 0 < a < 3/2, which covers 
a wide range of possible inner profiles. The cuspy profiles 
of a > 0 are motivated by the results of IV-body simula¬ 
tions. The profile with a = 1 corresponds to the so-called 
universal profile proposed by NFW (see also Lokas & Ma- 
mon 2001 ) as a fit to the profiles of simulated haloes, while 
the profile with a = 3/2 is identical to the one following 
from higher resolution simulations of Moore et al. (1998). 
The profile with a = 0, possessing an inner core is favoured 
by observations of spiral galaxies (see the Introduction) and 
is very similar (but not identical) to the profile proposed by 
Burkert (1995). 

The scale radius r s introduced in equation (|l7]) marks 
the distance from the centre of the object where the slope 
of the profile is the average of the inner and outer slope: 
r -( 3 +o )/2 Additional parameter which controls the shape 
of the profile is the concentration which we define here in a 
traditional way as 


c = 


(18) 


where r v is the virial radius, i.e. the distance from the centre 
of the halo within which the mean density is v = 200 times 
the present critical density, p C rit,o- The density of 200p C rit,o 
is traditionally assumed by most IV-body simulators to be 
the mean density of the object collapsing and virializing at 
the present epoch. This value, following from the spherical 
collapse model, depends however on the underlying cosmol¬ 
ogy, e.g. v = 178 is valid for the Einstein-de Sitter model, 
while we have v ~ 100 — 120 for the currently most pop¬ 
ular A(Q)CDM models with SIm = 0.3 and ^Ia(Q) = 0.7 
(Lokas & Hoffman 2001a, 2001b). However, since this study 
is motivated by the dark halo profiles obtained in A-body 
simulations we keep v = 200 in order to have the same defi¬ 
nition of concentration parameter as in most simulations. 

We normalize the density profile ( |l7| ) so that the mass 
within r v is equal to the so-called virial mass 


4 o 

M v = -7rr 3 np C rit,o- (19) 

The characteristic density of equation (|l7|) then becomes 


„ _ (3 

pchar — 


a)upcrit,o c“ 

3 F(C) 


( 20 ) 


where F(c ) is given by the hypergeometric function 


F(x) — 2 Fi( 3 — a, 3 — a; 4 — a; — x). 


( 21 ) 


The dark mass distribution following from (|l7]), ( fl9| ) and 

© is 

M d ( S ) = M„ S 3 -“^M ( 22 ) 

where we introduced s = r/r v . 

The concentration of simulated dark matter haloes has 
been observed to depend on the virial mass. Jing & Suto 
(2000) tested the relation c(M v ) for the masses of the order 
of normal galaxies and clusters in the case of density pro¬ 
files with a = 1 and a = 3/2 and found the concentration 
to decrease with mass with the highest c~ 10 (a = 1 ) and 
c ss 5 (a = 3/2) observed for the mass scale of the order of 
a normal galaxy (10 12 h~ 1 Mq). In the previous work (Lokas 
2001) we used an extrapolation of the c(M v ) relation fitted 
to the simulations results to smaller masses characteristic 
of dwarf galaxies. Here, however, similar relation would be 
needed for the a = 0 profile and the information provided 
by observations is very scarce. Although relations e.g. be¬ 
tween the inner mass and core radius have been found in 
the case of Burkert profile (Burkert 1995) and those can be 
translated into an c(M„) relation useful here, the estimates 
of this kind have been performed only for a small number 
of spiral galaxies and the errors of the fitted parameters are 
large (Borriello & Salucci 2001). We have therefore decided 
to assume different concentration parameters c = const in 
the further analysis. If we believe, however, that the trend 
of c decreasing with M v extends to smaller masses (and the 
relation was found to be even steeper for subhalos by Bul¬ 
lock et al. 2001) we should expect the cuspy halos in dSpli 
galaxies to have c > 10 (a = 1) and c > 5 (a = 3/2). 


4 RESULTS FOR THE FORNAX AND DRACO 
DWARFS 

In the two following subsections we model the velocity dis¬ 
persion profile of Fornax and Draco first assuming isotropy 
and then for arbitrary anisotropy parameter (3= const. In the 
third subsection we present predictions for the line-of-sight 
kurtosis for some of the best-fitting models. 


4.1 Isotropic orbits 

The observational parameters of the Fornax and Draco 
dwarfs needed in the following computation were taken from 
Irwin & Hatzidimitriou (1995) and are summarized in Ta¬ 
ble []} The Table gives the fitted Sersic radius Rs assum¬ 
ing the m = 1 Sersic distribution of equation , the dis¬ 
tance modulus to — M, the distance in kpc with its error, 
the brightness of the galaxies in V-band, My, and the to¬ 
tal luminosity in V-band, Ltot.v, together with its error. We 
adopt the mass-to-light ratio for stars in this band to be 
Ty « 1 Mq/Lq (Mateo et al. 1991). 

In the previous work (Lokas 2001) it was shown that the 
density of stars alone cannot produce the observed velocity 
dispersion profiles in the case of Fornax and Draco. In the 
following we will explore the hypothesis that the velocity 
dispersions are generated by stars moving in the Newtonian 
gravitational acceleration generated by stars and dark mat¬ 
ter. 

We first assume that the velocity distribution of stars 
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Figure 1. Upper panel: the best fitting masses of dark haloes as 
a function of assumed concentration parameter for profiles with 
different a in the case of isotropic orbits for the Fornax dwarf. 
Lower panel: the goodness of fit measure, % 2 /N , for the fits in 
the upper panel. 


Table 1. Observational parameters of the Fornax and Draco 
dwarfs from Irwin &; Hatzidimitriou (1995). 


parameter 

Fornax 

Draco 

Rs (aremin) 

9.9 

4.5 

Rs(kpc) 

0.35 

0.094 

m — M 

20.4 

19.3 

distance (kpc) 

120 ±8 

72 ±3 

My (mag) 

-13.0 

-8.3 

Ltot, v(Lq) 

(1.4 ±0.4) x 10 7 

(1.8 ±0.8) X 10 5 




Figure 2. Same as Figure 111, but for Draco. 


is isotropic, i.e. we consider models obtained from equations 
(^) and (^) in the case of [3 = 0. The 3D distribution of 
stars is(r) is given by (§)■ The gravitational acceleration is 
g = —d3>/dr = —GM(r)/r 2 with M(r ) = M*(r) ± 
where M*(r) and Md(V) are given by equations (^6|) and 
(p2[), respectively. Using such modelling of the velocity dis¬ 
persion profile we perform a one-parameter fitting of the ve¬ 
locity dispersion data. In the case of Fornax we use the data 
of Mateo (1997) and in the case of Draco those of Arman- 
droff, Pryor & Olszewski 1997 (based on the observations 
of Armandroff, Olszewski & Pryor 1995). In each case we 
assume the value of the concentration parameter c and fit 
by the method of least squares (Brandt 1999) the dark virial 
mass M v of the models with different dark matter profiles 
distinguished by the parameter a of equations (|l7j) and FL 
The resulting best-fitting values of M v as a function of 
assumed c are shown in the upper panels of Figures |l] and 
H The error bars indicate 1 <t errors in M v associated with 
the errors in the velocity dispersion measurements. They 
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a=3/2 isotropic orbits 
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Figure 3. Upper panel: the best fitting velocity dispersion pro¬ 
files obtained for c = 30 (« = 0), c = 20 (a = 1) and c = 10 
(a = 3/2) in the case of isotropic orbits for the Fornax dwarf. 
Lower panel: mass-to-light ratios for the fits in the upper panel. 



Figure 4. Same as Figure but for Draco. The best fitting 
velocity dispersion profiles plotted were obtained for c = 50 (a = 
0), c = 20 (a = 1) and c-- 10 (a - 3/2). 


are obtained from the covariance matrix (in this case one- 
element) of the fitted parameters. The errors in mass in the 
whole range of concentrations shown in Figures [l] and are 
of the order of 16-45% for Fornax and 24-185% for Draco. 
The lower panels of these Figures plot /N, the measure 
of the goodness of fits presented in the upper panels, where 
N is the number of degrees of freedom (in this case we have 
TV = 13 for Fornax and IV = 4 for Draco). The results for 
a = 0 do not reach to very low concentration because the 
least-square fitting converges more and more slowly there 
yielding incredibly high dark masses. One might think that a 
two-parameter (c, M v ) fitting would be a better choice here, 
but as the lower panels of Figures [I] and ^ show, this would 
give the lowest possible c in all cases. In the case of cuspy 
profiles resulting from 7V-body simulations low c of the order 
of unity is never obtained except for objects of the mass 
of cluster of galaxies. Clearly for both Fornax and Draco 


the a = 0 profile provides the best fit in terms of /N. 
Although quite low y? /N can be also obtained for a = 1 
profile, this happens only for unrealistic low concentrations 
and yielding incredibly high mass of a halo of the order of a 
normal galaxy. 

The predictions for the line-of-sight velocity dispersion 
profiles of the best-fitting models with a = 0,1 and 3/2 in 
the isotropic case are shown in the upper panels of Figures || 
and for Fornax and Draco respectively together with the 
data. For clarity of the illustration we have chosen only one 
of the fits from Figures [I] and |2j for each type of profile. In 
the case of a = 0 profile no a priori information on c is 
available so we have chosen the value giving the best fit and 
reasonably low best-fitting mass M v at the same time, i.e. 
c = 30 in the case of Fornax and c = 50 in the case of Draco. 
For cuspy profiles the concentrations expected on the basis 
of the IV-body results were chosen (see Subsection 3.2). We 
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Table 2. Examples of best-fitting parameters in the isotropic 
and anisotropic case with 1 <t error bars. M v is given in units of 
10 9 M q . 


case 

a 

parameter 

Fornax 

Draco 

(3 = 0 

0 

M v 

3.8 ± 1.0(1.2) 

4.1 ± 3 . 0 ( 3 . 2 ) 


1 

M v 

1.5 ± 0.3(0.3) 

1.9 ± 0.8(0.8) 


3/2 

M v 

1.3 ± 0.3(0.3) 

0.84 ±0.27(0.27) 

0# 0 

0 

M v 

3.4 ± 1.0(1.1) 

4.3 ± 3.4(3.6) 



(3 

-0.36 ±0.33 

-0.11 ±0.47 


1 

My 

1.3 ± 0.3(0.3) 

2.4 ± 1.1(1.2) 



(3 

-1.5 ±0.8 

-0.98 ± 1.2 


3/2 

M v 

1.2 ± 0.3(0.3) 

1.2 ± 0.4(0.4) 



(3 

-2.6 ± 1.6 

-3.5 ± 5.1 


have taken c = 20 for a = 1 and c = 10 for a = 3/2. The 
best-fitting virial masses obtained for these concentrations 
in the isotropic case together with lcr error bars following 
from the uncertainties in velocity dispersion measurements 
are given in the upper part of Table The method of least 
squares applied here allows also to study the propagation of 
errors associated with other observables entering the mod¬ 
els. Including the errors in luminosity and distance given in 
Table [j] and assuming an error in the stellar mass-to-light ra¬ 
tio of 0.5 Mq/Lq (an estimate given by Mateo et al. (1991) 
in the case of Fornax) we obtain errors for M v which are 
given in brackets in Table 

Lower panels of Figures |5| and ^ show the combined 
mass-to-light ratio in V-band for Fornax and Draco calcu¬ 
lated from 


M/Ly 


M D (r)+M.(r) 

Ly{r) 


(23) 


where Mrt(r) and A7*(r) are given by equations J 22 I ) and 
© respectively, and Ly(r) = M*(r)/ Ty is the luminosity 
distribution. Calculation of M/Ly were made for the same 
models as presented in the upper panels of the Figures. The 
growth of M/Ly at large distances is due to the behaviour 
of Mp(r) which diverges logarithmically for all considered 
a while the mass in stars is finite. However, while the cuspy 
profiles produce divergent M/Ly at the centre of the galaxy, 
the a = 0 profile gives a finite M/Ly increasing at all radii. 

All the mass-to-light estimates shown give a similar 
M/Ly « 10 Mq/Lq in the case of Fornax and M/Ly « 
IOOA/q / Lq in the case of Draco at a scale of the order of 
two Sersic radii and are consistent with the values of the 
central mass-to-light ratios given by Mateo et al. (1991) and 
Armandroff et al. (1997). 

Kleyna et al. (2001) estimated mass-to-light ratio for 
Draco using their new velocity measurements (Kleyna et al. 
2002). Assuming isotropic orbits and the velocity dispersion 
to be constant or a linear function of distance they calcu¬ 
lated the mass needed to produce such velocities from the 
Jeans equation. The resulting M/Ly at three core radii (30 
arcmin) turned out to be in the range of 350-1000. For this 
distance our isotropic models give M/Ly = 593, 369 and 250 
for a = 0,1 and 3/2 respectively. 


4.2 Anisotropic orbits 

Comparing the predictions shown in upper panels of Fig¬ 
ures |i] and ^ with the data we find that the cuspy pro¬ 
files have to be aided by a certain amount of tangential 
anisotropy in order to better reproduce the shape of the 
observed velocity dispersion profiles. Trying more radial 
anisotropy than in the case of isotropic orbits by using posi¬ 
tive values of (3 makes the curves decrease even more steeply 
with distance, contrary to the trend observed in the data. 
In the case of the profile with a = 0 it is not clear whether 
more radial or more tangential anisotropy would provide 
a better fit. Therefore in what follows we perform a two- 
parameter fitting of the velocity dispersion profiles with ve¬ 
locity anisotropy (3 =const as a second free parameter in 
addition to the dark mass, M v . 

The results of such least squares fitting procedure of the 
Fornax and Draco data are shown in Figures || and The 
Figures are analogous to Figures [I] and [ 2 ] for the isotropic 
case except for the added uppermost panels showing the 
best-fitting (3 parameters with their lcr error bars. The es¬ 
timated dark virial masses in the middle panels are now 
roughly the same for Fornax and somewhat higher for Draco 
than for isotropic orbits while their errors are now 17-50% 
for Fornax and 31-105% for Draco for the whole range of 
concentrations. The quality of fits (bottom panels) in terms 
of X 1 /N (here N = 12 for Fornax and N = 3 for Draco) is 
now significantly better for cuspy profiles and also somewhat 
improved for the profiles with a = 0. However, the profiles 
with cores produce best-fitting (3 closest to zero. 

Figures R and are the anisotropic analogues of Fig¬ 
ures and ^ As in the isotropic case for a = 0 profile we 
have chosen c = 30 in the case of Fornax and c = 50 in the 
case of Draco, while for cuspy profiles we have taken c = 20 
for a = 1 and c = 10 for a = 3/2. These values of concentra¬ 
tions happen to be those where yf /N reaches minimum for 
different profiles in the case of Fornax (see the lowest panel 
of Figure |). Upper panels of Figures [?] and ^ prove that 
when arbitrary anisotropy parameters are allowed accept¬ 
able fits can be obtained with all profiles. However, cuspy 
profiles require significantly more tangential anisotropy, i.e. 
lower (3 than the profile possessing a core. 

The best-fitting values of M v and (3 used in Figures ^ 
and H are given in the lower part of Table [ 2 ] together with 
their errors following from the uncertainties in velocity dis¬ 
persion measurements. As in the isotropic case the errors 
in brackets for M v include the effect of uncertainties in lu¬ 
minosity, distance and stellar mass-to-light ratio. The error 
estimates of (3 are not significantly affected by these addi¬ 
tional factors and are not listed. It is clear from Table ^ 
that in all cases the measurements of velocity dispersion are 
the dominant source of error, while other sources affect the 
results only slightly and more for dark matter profiles with 
lower a. One may object, however, that the distance errors 
given in Table ill are rather small and do not encompass 
higher distance estimates by other authors (Mateo 1998b). 
Increasing the adopted distances bv 10 kpc would decrease 
the virial masses given in Table El by up to 10% for For¬ 
nax and 40% for Draco for a = 0 with almost no difference 
for a = 3/2 (the estimates of (3 would remain roughly the 
same). 

One may also wonder how the best-fitting parameters 
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Figure 5. The best fitting anisotropy parameters (3 (upper panel) Figure 6. Same as Figure but for Draco. 

and masses of dark haloes (middle panel) as a function of assumed 

concentration parameter for profiles with different a in the case of 

anisotropic orbits for the Fornax dwarf. Lower panel: the goodness 

of fit measure, x 2 > f° r the fits in the upper panels. 
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Figure 7. Upper panel: the best fitting velocity dispersion pro¬ 
files obtained for c = 30 (a = 0), c = 20 (a = 1) and c = 10 
(a = 3/2) in the case of anisotropic orbits for the Fornax dwarf. 
Lower panel: mass-to-light ratios for the fits in the upper panel. 


are affected by our definition of the virial mass. It turns out 
that decreasing the parameter v in equation (juj) by a factor 
of two changes mainly best-fitting M v (leaving j3 roughly 
the same) for the same value of concentration by increasing 
it typically by a factor of a few, more for lower a and lower 
c (it must be emphasized however, that c given by equation 
( [fq ) would then be a different quantity since r v for a given 
M v would be increased). 

The mass-to-light ratios shown in the lower panels of 
Figures [t| and ^ do not change significantly with respect to 
the isotropic case and the very different behaviour for cuspy 
profiles versus the profile with a core is preserved. In the 
case of Draco at R = 30 arcmin our anisotropic models give 
M/Ly = 606,414 and 300 for a = 0,1 and 3/2 respectively, 
values slightly higher than in the isotropic case. The values 
agree to the order of magnitude with the estimate given by 




Figure 8. Same as Figure [J, but for Draco. The best fitting 
velocity dispersion profiles plotted were obtained for c = 50 (a = 
0), c = 20 (a = 1) and cm 10 (a = 3/2). 


Kleyna et al. (2002) who find mass-to-light ratio of 330 for 
their best-fitting anisotropic model of Draco. 


4.3 Predictions for the kurtosis 

In the previous subsection we have shown that when velocity 
anisotropy is treated as a free parameter equally good fits to 
velocity dispersion profiles can be obtained with both cuspy 
profiles and the profile with a core. In order to get rid of this 
degeneracy and discriminate between the profiles one has to 
resort to higher order moments of the velocity distribution 
such as the kurtosis introduced in Section 2. Figure^ shows 
the predictions for this quantity for models which turned out 
to provide best fits to velocity dispersion profiles of Fornax 
and Draco. The parameters used for every type of profile 
were taken from Table In each panel the thick lines show 
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Figure 9. Kurtosis predictions for the best fitting models of Ta¬ 
ble |2j for Fornax (upper panel) and Draco (lower panel). Thick 
lines show the results for anisotropic models while thin lines for 
isotropic ones. 


the predictions for anisotropic models while thin lines for 
isotropic ones. 

The results prove that the kurtosis can indeed differen¬ 
tiate between the profiles with a cusp and those with a core. 
In the anisotropic case we find for both galaxies that the 
profile with a core results in a higher value of the kurtosis 
at a given distance from the centre and its value decreases 
more steeply with R than in the case of cuspy profiles. The 
situation is similar in the case of isotropic orbits, but here 
«ios(-R) definitely decreases with distance only for the a = 0 
profile while its values remain almost constant for the cuspy 
profiles (slightly decreasing for a = 1 and slightly increasing 
for a = 3/2). 


Table 3. Best-fitting parameters for MOND in the isotropic 
and anisotropic case with I a error bars, ao is given in units of 
10“ 8 cm/s 2 , while M/Ly in Mq/Lq. 


case 

assumed 

fitted 

Fornax 

Draco 

d = 0 

M/Ly = 1 

a 0 

2.2 ± 0.5(0.8) 

39 ± 13(22) 


ao = 1.2 

M/Ly 

1.8 ± 0.4(0.7) 

32 ± 10(18) 


M/Ly = 1 

ao 

2.1 ± 0.5(0.8) 

50 ± 19(29) 



a 

-0.96 ±0.53 

— 1.5 ± 1.6 


ao = 1.2 

M/Ly 

1.8 ± 0.4(0.7) 

40 ± 14(23) 



0 

-0.97 ±0.53 

— 1.5 ± 1.7 


5 ALTERNATIVE: MOND 

DSph galaxies have low densities and therefore they lie in 
the regime of small accelerations. According to Modified 
Newtonian Dynamics (MOND), the theory proposed by Mil- 
grom (1983), for accelerations below the characteristic scale 
of ao s» lCP 8 cm/s 2 the laws of Newtonian dynamics may 
change affecting either inertia or gravity and mimic the pres¬ 
ence of dark matter. Although originally designed to explain 
flat rotation curves of spiral galaxies, MOND is expected 
to operate in all low-acceleration systems. It has been e.g. 
shown to explain many observed features of LSB galaxies by 
McGaugh & de Blok 1998. On the other hand, Lake (1989) 
claimed MOND to be inconsistent with rotation curves of a 
few dwarf galaxies (the required values of a o being different 
for different systems and generally too low). In the case of 
dSph galaxies MOND should cause an increase of velocity 
dispersions. 

It is presently far from clear, whether MOND can ex¬ 
plain the observed velocity dispersions of dSph galaxies. Al¬ 
though for most dwarfs it gives satisfactory results, noto¬ 
rious problems have been caused by Draco. Analyzing its 
central velocity dispersion Gerhard & Spergel (1992) have 
found it incompatible with stellar mass-to-light ratios, con¬ 
trary to the conclusion of Milgrom (1995) who claimed that 
its mass-to-light ratio is consistent with stellar values if all 
the observational errors are properly taken into account. 

In the previous work (Lokas 2001) we have shown how 
to generalize the method of predicting the line-of-sight ve¬ 
locity dispersion profile to MOND and applied it to fit the 
data for Fornax, Draco and Ursa Minor dwarfs. We have as¬ 
sumed the stellar mass-to-light ratio of M/Ly = IMq/Lq 
and fitted parameters ao and f3. The best-fitting values of 
ao were different for different dwarfs, acceptable for Fornax, 
but much higher than expected in the case of Draco and 
Ursa Minor. We concluded that while the case of Ursa is 
doubtful since it has been recently observed to possess tidal 
tails (Martfnez-Delgado et al. 2001), the high value of ao 
estimated for Draco can pose a serious problem for MOND 
(a similar conclusion based on the analysis of a velocity dis¬ 
persion profile of Draco from a different data set was also 
recently reached by Kleyna et al. 2001). 

We have repeated the analysis here for Fornax and 
Draco with the least-square fitting method, similar to that 
applied to fit dark matter profiles, in order to provide the 
error estimates for the fitted parameters. Both isotropic and 
anisotropic cases were studied. The results of the fitting pro¬ 
cedure are shown in Table with the isotropic case in the 
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upper part and the anisotropic in the lower. We assumed 
M/Ly = 1 Mq/Lq and fitted the acceleration parameter 
oo in one case while taking ao = 1.2 x 10 ~ 8 cm/s 2 , the value 
known to fit best the rotation curves of spiral galaxies (Bege- 
man, Broeils & Sanders 1991; Sanders & Verheijen 1998), 
and finding M/Ly in the other. 

As is clear form Table ^ both approaches yield accept¬ 
able values of M/Ly or ao in the case of Fornax, while very 
high numbers in the case of Draco. However, the errors for 
Draco estimated from the uncertainties in the velocity dis¬ 
persion measurements are so large that the acceptable val¬ 
ues of M/Ly or ao (of the order of unity in the units of 
the Table) are still within 3<r error bars. Analysis of the 
velocity dispersion profile therefore reduces the error sig¬ 
nificantly compared to what can be obtained from a sin¬ 
gle central velocity dispersion measurement (then accept¬ 
able M/Ly are within 2<r error from the best-fitting value, 
McGaugh, private communication), but still not enough to 
rule out MOND. In fact these errors are farther increased 
if the errors in luminosity measurement (quite big in the 
case of Draco, see Table hi) and in distance are taken into 
account. The resulting larger errors (due mainly to the error 
in luminosity) are given in brackets in Table 111 With these 
errors the values obtained here are safely within 2cr of the 
acceptable ones. (A discussed in Section 4 this source of er¬ 
ror is not important for dark matter because it dominates 
strongly the distribution of stars.) We therefore conclude 
that with present quality of the data the alternative expla¬ 
nation of velocity dispersions in terms of MOND cannot yet 
be ruled out. 


6 DISCUSSION 

We presented predictions for the line-of-sight velocity disper¬ 
sion profiles of dwarf spheroidal galaxies and compared them 
to observations in the case of the Fornax and Draco dwarfs. 
We discussed different dark matter distributions, both cuspy 
and possessing flat density cores. For isotropic models the 
dark haloes with cores are found to generally fit the data bet¬ 
ter than those with cusps. The fits for each type of profile 
tend to be better for lower concentrations. For cuspy profiles 
those cannot be taken arbitrarily low since, as discussed in 
Subsection 3.2, N-body simulations (providing the basis for 
the choice of profiles here) predict higher concentrations for 
smaller masses. Besides, fits with low concentration require 
incredibly high corresponding dark masses of the order of a 
normal galaxy. In our exemplary cases we have taken c = 10 
for q = 3/2 and c = 20 for a = 1. The dark matter profiles 
with such parameters yield for Draco fits with y/ 2 /N = 1.6 
and \ 2 /N = 1, respectively, while an a = 0 profile with 
c = 50 has yf /N = 0.6. The situation is similar for Fornax. 

Anisotropic models were studied by fitting two param¬ 
eters, dark mass and velocity anisotropy, to the data. We 
find that the steeper the cusp of the profile, the more tan¬ 
gential is the velocity distribution required to fit the data, in 
agreement with the well-known degeneracy of density pro¬ 
file versus velocity anisotropy. In this case acceptable fits 
(with x 2 /N = 1 or lower) can be obtained for all considered 
profiles, but only the profiles with cores yield good fits with 
anisotropy parameter close to the isotropic value. To dis¬ 
criminate between different profiles we offer predictions for 


the kurtosis of the line-of-sight velocity distribution. Our 
findings in the case of the velocity dispersion profile for 
Draco are complementary to and in qualitative agreement 
with those of Kleyna et al. (2002) who tested dark matter 
distributions varying at large distances. Their best-fitting 
models required the halo to be nearly isothermal (as our 
models are in a wide range of distances) and the velocity 
anisotropy to be weakly tangential. 

The amount of dark matter inferred depends on the 
concentration parameter of the density profile, but for real¬ 
istic concentrations is typically of the order of 10 9 M,V) both 
for Fornax and Draco for all profiles and velocity distribu¬ 
tions. Since the total luminosity of Draco is two orders of 
magnitude lower than that of Fornax (see Table ^), this re¬ 
sults in a significant difference in their corresponding mass- 
to-light ratios. At the scale of the order of two Sersic radii 
Fornax has M/Ly as 10 Mq/Lq while for Draco we obtain 
M/Ly as 100 Mq/Lq. However, mass-to-light ratio depen¬ 
dence on distance is very different for different profiles. For 
cuspy profiles M/Ly diverges both at small and at large dis¬ 
tances from the centre, while for a = 0 profile it is finite at 
the centre and grows steadily with distance. For distances 
larger than two Sersic radii profiles with lower a give higher 
values of M/Ly and the anisotropic models yield slightly 
higher M/Ly than isotropic ones. As discussed in Section 
3, our results for M/Ly outside the galactic core agree to 
the order of magnitude with the estimates of Kleyna et al. 
( 2001 , 2002 ). 

A possible source of uncertainty in the analysis of veloc¬ 
ity dispersions is the influence of binary population. In the 
recent study De Rijcke & Dejonghe (2002) using realistic 
distributions of orbital parameters estimated the maximum 
additional velocity dispersion due to binary stars to be of 
the order of 3 km/s (in agreement with earlier studies). It 
turns out that only stellar systems with intrinsic dispersions 
close to this value (e.g. globular clusters) will be affected. For 
velocity dispersions characteristic of dSph galaxies studied 
here the effect is small: the velocity dispersion a = 9 km/s 
for a galaxy with a Gaussian line-of-sight velocity distribu¬ 
tion can increase due to the presence of binaries by at most 
5% (for the maximum binary fraction of 1, see Figure 15 of 
De Rijcke & Dejonghe 2002), much less than the typical er¬ 
ror in velocity dispersion measurements. The binary fraction 
of Draco is estimated to be between 0.2 (Olszewski, Pryor & 
Armandroff 1996) and 0.4 (Kleyna et al. 2002), which would 
make the effect even smaller. For Fornax no information on 
binaries is yet available, but the typical velocity dispersions 
are higher in this galaxy which makes the effect less im¬ 
portant. We therefore expect binaries to only slightly affect 
our best-fitting models so as to increase the estimated virial 
mass. More influence is expected in the case of kurtosis: its 
Gaussian value Ki os = 3 for the same dispersion a = 9 km/s 
will increase by up to 23%. Hence when using the kurtosis 
predictions of Subsection 4.3 to discriminate between differ¬ 
ent dark matter profiles the influence of binary population 
will have to be studied in more detail. 

The analysis presented in this paper is valid under the 
condition that dSph galaxies are in dynamical equilibrium. 
However, those of them studied here, as members of the 
Local Group exist in gravitational field of larger galaxy, 
the Milky Way. One may worry that the tidal interactions 
with bigger galaxy may affect the dynamics of some dwarfs 
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and our interpretations of their velocity dispersions. Some of 
the dwarfs have been indeed observed to possess tidal tails, 
e.g. Ursa Minor (Martfnez-Delgado et al. 2001), Sagittarius 
(Ibata, Gilmore & Irwin 1994) or Carina (Majewski et al. 
2000). However, in the case of Draco dwarf studied here no 
evidence of tidal interactions was found (Piatek et al. 2002). 

Tidal interactions have been studied theoretically by 
numerical modelling with conflicting results. Piatek & Pryor 
(1995) have performed numerical simulations of such effects 
and concluded that they should not affect much the inferred 
mass-to-light ratios. However, Kroupa (1997) and Klessen & 
Kroupa (1998) have demonstrated that what appears to be 
a dSph galaxy in equilibrium may in fact be a tidal remnant 
displaying high velocity dispersions without any dark mat¬ 
ter. As suggested by Oh, Lin & Aarseth (1995) it may also 
very well be that some dwarfs indeed contain dark matter 
while others may be part of tidal debris. It is still debated, 
whether tidal interactions are common among dSph galax¬ 
ies and whether they may cause an increase of velocity dis¬ 
persion that is likely to be interpreted as presence of dark 
matter. 
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